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We present results of molecular dynamics simulations on lithium metasilicate over a broad 
! ■ range of temperatures for which the silicate network is frozen in but the lithium ions can still be 

■ equilibrated. The lithium dynamics is studied via the analysis of different correlation functions. 
, The activation energy for the lithium mobility agrees very well with experimental data. The 

l/-) \ correlation of the dynamics of adjacent ions is weak. At low temperatures the dynamics can 
^ | be separated into local vibrational dynamics and hopping events between adjacent lithium 

■ sites. The derivative of the mean square displacement displays several characteristic time 
O ■ regimes. They can be directly mapped onto respective frequency regimes for the conductivity. 
"£3 . In particular it is possible to identify time regimes dominated by localized dynamics and long- 
£h ! range dynamics, respectively. The question of time-temperature superposition is discussed for 
^ \ the mean square displacement and the incoherent scattering function. 

: 

o . 

P. 1 Introduction 



^ , The dynamics of ions in amorphous materials is very complex, as indicated, e.g., by the strong 
frequency dependence of the conductivity [0], |j. This frequency dependence directly 

reflects the presence of back- and forthdynamics of the individual ions. The complexity of the 
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ion dynamics is due to the simultaneous action of the time- dependent Coulomb interaction 
with the other mobile ions and the time-independent interaction with the spatially disordered 
and basically immobile network. Interestingly, the frequency dependence of the conductivity 
is very similar when comparing different materials M. Even simulations of disordered hopping 
models display a similar frequency dependence; see, e.g., |5|. |6|. 

Typically the frequency-dependence of the conductivity a{y) displays several characteristic time 
regimes. For a sodium silicate system |7j one observes, e.g., a microsopic regime for v > 10 13 
Hz, a scaling regime <r(z/) oc v 2 for 10 13 Hz > v > 10 11 Hz, and finally for v < 10 11 Hz 
a continuous decrease of cx(z/) with decreasing v. In a double-logarithmic representation at 
low temperatures the apparent exponent decreases from one to zero until the d.c. plateau is 
reached, i.e. <y{y) = Od.c- A possible interpretation of this behavior relates these frequency 
regimes to local vibrations, stochastic localized dynamics, and jump dynamics, respectively 
|§. An important observation for the frequency-dependent conductivity uiv) is the validity of 
the time-temperature superposition principle. For most materials it can be expressed by the 
Summerfield scaling for which a{p)/ad c is only a function of vjiTodcj li- 
lt is important to relate the frequency dependence of a{y) to the real space dynamics of the 
ions. In this way one may get interesting information about the nature of the ion dynamics. 
This can be done by expressing in terms of the mean square displacement (r 2 (t)). Con- 
ceptually, both observables differ in one important aspect. While the conductivity contains 
effects of multi-particle correlations, the mean square displacement is a single-particle quantity. 
Experimentally, the correlation is expressed by the Haven ratio []T0] , |TT| which in general can be 
frequency dependent (denoted Hr{v)). The experiments indicate that the correlation among 
different ions is small and furthermore can be approximated by a value Hr. 1/Hr is a measure 
for the number of ions which are significantly correlated. Then one can relate er(z/) to (r 2 (t)) 
via linear response theory [[[], Q 

a 2 f°° 

a{v) = — H dt(d/dt)w{t)cos{2irvt) (1) 

bH R pk B -L Jo 

where q denotes the charge, p the density of the mobile ions, and Hr the Haven ratio. The 
function w(t) is defined as 

w(t) = (d/dt)(r 2 (t)). (2) 

Analysis of spatial aspects of dynamical processes is one of the strongholds of computer sim- 
ulations. Here it is of particular help that detailed information is available on a microscopic 
level. Conceptually, one can proceed in two steps. In a first step one can study w(t). On the 
basis of Eq.[l] it has been readily shown that for a significant dispersion of the conductivity 
the frequency-dependence of cr(i>) and the time-dependence of w(t) are basically mirror images 
of each other when identifying t with \j-K 2 v [0]. Therefore it is the quantity w(t) which is of 
uttermost importance for the direct relation to conductivity experiments since one may hope 
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to recover the three regimes of cr(y) also in w(t). As a consequence of Eq.l a dispersion in 
cr(v) is equivalent to a time dependence of w(t) and a scaling relation w(t) oc t~ a transfers 
to cr(u) oc u a . In a second step one can analyse the dynamics in more detail to elucidate the 
underlying physical processes in the different time regimes. Although interesting details may be 
seen from the analysis of individual trajectories a sound understanding of the dynamics should 
be based on the study of appropriate correlation functions, yielding the average dynamical 
behavior. 

In this paper we present computer simulations for lithium silicate following this general ap- 
proach. As an example we take the lithium dynamics in (Li20)(Si02) |T3|, H, 13]. In contrast 
to most earlier work on this system the present computer generation allows one to study the 
lithium dynamics in the low-temperature dispersive regime in equilibrium conditions. Sim- 
ilar systems like sodium silicate systems are also of current interest for numerical studies 
|l6i |i~7l , |i~8] , [19| |20| , |2l|, p2j [23t [24], [25j. The outline of this paper is as follows. In Section 
2 we describe the technical aspects of the simulation and discuss the numerical tools. The 
results of our simulations as well as their interpretation are presented in Section 3. We close 
with a discussion and a summary in Section 4. 



2 Simulation 

The potential energy of the lithium silicate system can be written as the sum of a Buckingham 
and Coulomb pair potential 

Unix) = ^jr- ~ 7? + A » «p(-%r). (3) 

The pair indices i,j characterize the various ion pairs (Li-Li, Li-Si, etc.). The individual potential 
parameters are listed in our previous work [E(J . They are based on the work of Habasaki et al. 



13], |14| , p27| , Pq , [29f] and have been obtained from ab-initio calculations. Effective charges qsi = 
2.4, qu = 0.87, and qo = 1-38 have been chosen. This choice fulfills charge neutrality of the 
total system. Periodic boundary conditions have been used. Simulations during the production 
runs were performed in the NVE ensemble. The elementary timestep of our molecular dynamics 
simulations was 2 fs, the density p = 2.34 g cm" 3 , taken from experimental room temperature 



data PU[ . This corresponds to pressures of the order of 1 GPa at the lowest temperatures of our 



simulation. The total number of atoms in our system is 1152. For generating the trajectories we 



have used the programm MOLDY, supplied by K. Refson f31fl . In our longest simulation run the 
system has been propagated for 20 ns. For our simulations we have started from configurations 
at T = 1500 K where ions as well as network could be equilibrated. Equilibration of the total 
system at this temperature requires a simulation run of ca. 4 ns. We note that the computer 
glass transition as observed on the time scale of 10 ns is approximately 1100 K (see network 
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diffusion data in f26fl). In the present simulation we have chosen one temperature slightly 
above T g (T = 1240 K) and four temperatures below T g (T = 980 K, 750 K, 700 K, and 640 
K). For these four temperatures we have paid attention to equilibrate the lithium subsystem 
before recording our observables. The network only performs local fluctuations with an average 
radius much smaller than one Angstrom. Interestingly, despite this fact we observed a slight 
temperature drift as a result of the continuous relaxation of the network for temperatures below 
the glass transition. 



3 Results 

3.1 Trajectories 

In Fig.l some representative trajectories of individual lithium ions at T = 750 K are shown. 
For these ions the dynamics can be described as a series of local vibrations and jumps between 
adjacent ionic sites. Similar features have been observed in several previous simulations on 



alkali silicates, see e.g. Refs. |22|, |23], j24[ . The three ions, displayed in Fig.l, behave somewhat 
differently. Whereas ion (1) does not leave its intitial site at all, ion (2) performs several back- 
and forthjumps between two adjacent sites, and ion (3) follows a random-walk like path between 
adjacent sites. The variety of different motional patterns directly shows that it is important to 
use appropriate statistical observables to identify the average type of behavior of the lithium 
dynamics. This is the subject of the remainder of the paper. 

3.2 Structure 

In Fig. 2 we display the partial Li-Li structure factor g(r) for several temperatures. One can 
clearly see that g{r) shows a sharp nearest neighbor peak at r nn »s 2.6 A. Interestingly, the width 
of the nearest-neighbor peak of g(r) only slightly changes in the temperature range between 
T = 750 K and T = 1240 K. Here we have included also a very high temperature T = 4000 K 



from our previous simulations [3B|. Only for the higher temperature the nearest neighbor peak 
broadens. This shows that the low-temperature width of the nearest-neighbor peak is mainly 
due to the disorder of the network and not due to thermal fluctuations in crystalline-like ionic 
sites. 

To study the degree of clustering of the lithium ions we have calculated the function 

= Kdsg(s)4*s* 
y 1 47rr 3 /3 V ; 

h(r) counts the number of particles around a central particle in a sphere of radius r relative to 
the number one would have for a random distribution of particles. Thus the value of h(r) is a 
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direct measure whether or not a particle is surrounded by more (h(r) > 1) or less (h(r) < 1) 
particles than expected for a random distribution. The first neighbor shell, as obtained from 
the first minimum of g(r) ranges up to r « 4.4 A. Therefore the value h(r = 4.4 A) yields 
information whether the number of particles in the nearest neighbor shell is smaller or larger 
than expected for a statistical distribution. The function h(r) is displayed in Fig. 3. It turns out 
that for low temperatures h(r = 4.4 A) « 0.9. This indicates that the repulsion effect of the ions 
due to their mutual Coulomb interaction dominates the structure of the lithium subsystem. For 
T = 4000 K the distribution is essentially statistical, as expected for a dominance of entropic 
contributions. It is known experimentally from the study of alkali silicates |33, that for 
smaller alkali concentration the alkali ions start to form clusters and indeed we observe in so 
far unpublished work that for (Li 2 0)9(Si0 2 ) h(r = 4.4 A) « 1.5. 



3.3 Dynamics 

At first we would like to present results for the mean square displacement (r 2 (t)) at different 
temperatures; see Fig. 4. For all temperatures there exists a significant subdiffusive regime 
which is more pronounced for the lower temperatures. Note that at much higher temperatures 
the subdiffusive regime disappears and one observes a direct transition from a ballistic short- 
time regime (r 2 (t)) oc t 2 to a diffusive long-time regime (r 2 (t)) oc t |p6| . The diffusion constant 
D can be directly estimated from the long-time regime via D = (r 2 (t))/(Qt). 
Next we present the temperature dependence of the diffusion constant D(T). The results are 
shown in Fig. 5. We have included high-temperature data points as already presented in our 
previous work | 26|| . The dynamics is characterized by a single activation energy E a = 0.58 eV 



in the temperature range between T = 6000 K and T = 640 K. The computer glass transition 
T g « 1100 K does not show up in the temperature dependence of the diffusion constant because 
already above this temperature the ion and network dynamics are strongly decoupled. 
For a comparison of the activation energy with experimental results one should consider con- 
stant pressure rather than constant density simulations. In order to check how much the above- 
mentioned value of E a is modified we took the pressure of the T=1500 K run and checked how 
much the diffusion constant at T=6000 K is modified when using constant pressure rather 
than constant density. It turned out that the diffusion constant increased by a factor of 1.6. 
This corresponds to an increase of the activation energy by 0.08 eV, i.e. E a = 0.66 eV. This 
value compares very well with the experimental value of 0.59 eV as obtained from conductivity 



experiments [33 



In Fig. 6 we show that all curves can be superimposed on each other by substituting t by 
t* = (D(T)/D )t ( The constant D Q is chosen as D(T = 640 K)). Thus one may conclude that 
in the scaling regime, starting beyond the ballistic regime, (r 2 (t*)) fulfills the time-temperature 
superposition principle. The relation to the experimental findings will be discussed in the next 
Section. 
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From a theoretical perspective the incoherent scattering function S(q,t) is an important ob- 
servable to characterize the dynamics. It is defined as 



S(q,t) = (cos(q(f(t)-m)))- (5) 

For isotropic systems it only depends on the absolute value q of the wave vector q. In particular 
for q = q max = 2Ti/r nn (in our case q max = 2.4 A -1 ) one gets information about the dynamical 
processes on the lengthscale of nearest neighbor lithium distances. For glass forming liquids one 
can formulate a theory which characterizes the time and wave vector dependence of S(q, t) for 
a broad range of temperatures In Fig. 7 one can find S(q max ,t) at different temperatures. 
The decay time is an appropriate measure for the time scale on which an ion effectively moves 
to the next site. 

For glass-forming liquids it is predicted from theory and confirmed by simulation that S(q, t) 



follows the time-temperature superposition principle in the long-time regime |p5| , p6| . Inter- 
estingly, also for the present case of the lithium dynamics time-temperature superposition is 
observed. This can be seen from Fig.7 where at all temperatures S(q max ,t) is fitted to the 
Kohlrausch- Williams- Watts function Aexp(— (t/r) 13 ) with (3 = 0.45. Actually, the temperature 
dependence of r is also activated (not shown). In the range of temperatures analysed in this 
work the variation of r is four times larger than the variation of the diffusion constant. Thus 
translates into a higher activation energy of 0.70 eV rather than 0.58 eV as obtained for the 
diffusion constant. A possible origin of this descrepancy is discussed in the next Section. 
The primary goal of this paper is to elucidate the nature of the frequency dependence of the 
conductivity. As already discussed the conductivity and the mean square displacement are 
directly related if the Haven ratio Hr is constant. This is particularly the case if Hr is close 
to one, i.e. dynamical inter-ion correlations are small. In recent work we have introduced a 
quantitity N coop {t) which is a measure of cooperativity on a time scale t |fT7|] . It is defined as 



N v Y^xmm 

JWU - E . Xj(()Xj(() w 

where Xi(t) is a dynamic property of ion i on the time scale t which on average is zero, i.e. 
J2i.Xi(t) = 0. Here we choose 

x % {t) = {m-m)?-(r\t)). (7) 

For the practical implementation of this approach it is essential that N coop is calculated for 



subsystems. The technical details of this analysis can be found in |[37|| . The results are shown 
in Fig. 8. For a better comparison of the different temperatures we have plotted N coop against 
the normalized time t* . As expected N coop (t*) approaches one on very short and very long 
time scales since the different ions behave independently. For all temperatures the maximum 
occurs at approximately the same scaled time t* ~ 500ps. This time corresponds to the regime 
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where the transition from the strongly subdiffusive to the diffusive regime occurs. For the lower 
temperatures there are approximately 4 particles which are dynamically correlated. In contrast 
to glass-forming systems there is only a minor temperature dependence of N coop . A further 
difference for glass-forming liquids is that the maximum value of N coop can be larger than 20 
close to the glass transition j37]. Therefore the present values of N coop are indeed small. 



Is N coop related to the Haven ratio? On a qualitative level both quantities express the dynamic 
cooperativity. In the most simple case exactly M <C N particles behave identically but all 
these subsets of M particles behave in an uncorrelated manner. Then one simply has 

N coop = 1/H R = M. (8) 

On a more quantitative level we start from the zero-frequeny limit of the Haven-ratio, thereby 
slightly extending the arguments given in fl37| . This limit is given by 



R z itj j °° dt ( Vi (o) Vj (t)y {) 

These terms can be further rewritten in order to obtain the connection with N coop . Due to 
time-reversal symmetry of the molecular dynamics one has (vi(0)vj(t)) = (vi(0)vj(—t)). We 
define t m as a time far in the diffusive long-time regime, i.e. much longer as the time for 
which the velocity-correlations have basically decayed to zero. Then one can write 

/tm 
dt( Vi (0) Vj (t)). (10) 

For — t m +t d < t <t m — t d one can substitute ttj(0) by i>j(r). Thus, for t m — > oo one can write 

tm ftm ftm 

dt(v l (0)v J (t)) = (l/2t m )l dr dtiviir^it)). (11) 

-t m J tm i-m 

For stationary processes one can substitute the integration limits by and 2t m , respectively. 
Rewriting Eq.[|in terms of the above relations and replacing J 2 ' m dtvi(t) by rj(2t m ) — rj(0) one 
ends up by H R = l/N coop (t — > oo) when choosing Xi(t) = fi(t) — r*j(0). Thus Eq.|| also holds 
more generally. Unfortunately, for the present data the statistical error was very large with this 
choice of Xi(t) and much smaller with the choice in Eq.[7[ However, for the hard sphere system 
analysed in |J7| it turned out that to a very good approximation the maximum of N coop {t) 
with time for a quantity involving the absolute value of the shift vectors (like the one used 
above) is very close to the long-time limit of N coop for the choice Xi(t) = fi(t) —7^(0). Thus the 
maximum value of N coop in Fig. 8 is indeed a good measure for the inverse zero-frequency Haven 
ratio. Of course, a more detailed analysis of this aspect would be desirable. The resulting low- 



temperature estimate Hr fa 1/4 is in agreement with typical experimental values [11]. Finally 
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we would like to mention that even this modest cooperativity has been stressed as a key feature 
to understand the long-time diffusion [BS . 



For a closer discussion of w(t) (and thus of cr(u)) it will turn out to be helpful to use additional 
information about the real-space dynamics. For this purpose we have determined the van-Hove 
self-correlation function G s (r,t) for the lithium ions [B9| 



G s (r,t) = (S(r-\(r(t)-r(0))\))- (12) 

It denotes the probability that a lithium ion moves a distance r during the time t. The results 
are shown in Fig. 9. One can clearly see that at short times the dynamics of all ions is confined 
to a small r-range. Of course, for longer times the ions may explore larger regions of the system. 
The presence of the first well-resolved peak around do m 2.6 A shows that hopping processes 
with jump length do are a relevant feature of the dynamics. To a good approximation do is 
identical to the r-value r nn of the first nearest-neighbor peak of the partial structure factor 
g(r). Also the other peaks of the partial structure factor show up in G s (r,t). Thus one may 
conclude that the network provides well-defined lithium sites which are separated by barriers 
larger as compared to k^T . The lithium ions basically jump between adjacent sites. The term 
jump implies that the time scale it takes to cross the saddle is much shorter than the time scale 
the ions fluctuate within the individual sites. This conclusion is in qualitative agreement with 



previous simulations [15|, |24], Q and also agrees with the appearance of the trajectories 

in Fig.l. 

The first minimum of G s (r,t) around r min « 1.5 A can be used to distinguish local vibrational 
dynamics (|Ar(t)| < 1.5 A) and long-range dynamics (|Ar(t)| > 1.5 A). Of course, there may 
be individual ions which display local fluctuations with a length scale larger than 1.5 A. On 
average, however, this value of 1.5 A turns out as a significant length scale in the above- 
mentioned sense. This scenario of jumps between well-defined potential wells implies that the 
ions explore the individual sites for rather long times before they leave the site on a time scale 
approximately given by the decay time of S(q max , t). In order to analyse the intra-site dynamics 
somewhat closer we have determined the mean square displacement (r 2 (t))i oca i of those particles 
which at time t have moved less than 1.5 A. This curve is shown at all four temperatures of 
our study in Fig. 10. For times larger than 0.2 ps (r 2 (t))i oca i starts to become constant. This 
means that already on this short time scale the lithium ions have largely explored the local 
potential well. The increase of (r 2 (t))i oca i for T = 1240 K at long times is related to the fact 
that the network starts to relax since we are above its glass transition temperature. Note that 
the temperature dependence of this local exploration time is very weak. A careful inspection 
shows that for the two lower temperatures there is further gradual increase until t ~ 10 ps 
(T = 750 K) and t ~ 100 ps (T = 640 K), respectively. This observation indicates the presence 
of a substructure of the different wells. At low temperatures the exploration of the initial well 
might require the crossing of some lower saddles. This would explain the slower approach to 
the final plateau value as compared to the dynamics in a smooth local potential. 
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For long times all curves approach a constant plateau value (r 2 (t — ► oo))/ ocaZ . An exception is 
T = 1240 K where we take the value around t = 10 ps. In Fig. 11 this plateau value is plotted 
against temperature. To a good approximation (r 2 (t — > oo))i OC ai is proportional to temperature. 
This is expected for the dynamics in a harmonic potential. Thus the local dynamics can be 
characterized by an effective harmonic potential. These local fluctuations are related to the 
Debye- Waller factor ||. 

Now we are in a position to discuss w(t). The lengthscale r m i n = 1.5 A as derived from the 
van Hove self correlation function allows us to separate local and long-range motion and was 
the basis to study {rf occd (t)). In analogy to w(t) we define 

w local (t) = (d/dt){r 2 ocal (t)). (13) 

It is shown in Fig. 12 at different temperatures. The temperature dependence of wi oca i is rather 
weak. For t > 20 fs it roughly scales like t~ 2 . 

In Fig. 13 w (i) is shown for several temperatures. Three different time regimes can be distin- 
guished, (i) For t < 20 fs there is the ballistic regime, i.e. w(t) oc t. (ii) For 20/s < t < lps 
there is a strong decay of w(t). (iii) For t > lps w(t) decays much slower until it becomes 
constant for large times. Whereas for T = 1240 K and T = 980K this final decay is very weak 
it becomes significantly larger for lower temperatures. For T = 640 the decay spans even more 
than a decade. 

It is instructive to compare this with the typical experimental cr(i/) data for alkali silicates, as 
already discussed in the Introduction. The frequency regimes there have a one-to-one corre- 
spondence to the characteristic time regimes of w(t). (i) The vibrational modes for v > 10 13 Hz 
correspond to the ballistic regime in the time domain. Since l/(ir 2 v) ~ 10/s this agrees very 
well with the corresponding time regime of w(t) (ii) As mentioned before the <r{v) oc v 2 behav- 
ior for 10 13 Hz> v > 10 11 Hz can be attributed to stochastic dynamics in a harmonic potential 
||. For the time dependence of w(t) of an overdamped vibration one expects w(t) oc exp(— kt) 
with some decay constant k |[44|| . For a complicated system like an amorphous ion conductor a 
distribution of decay constants should be present, resulting in a somewhat modified decay char- 
acteristics, but still yielding a{y) oc v 2 for v < l/k max as obtained via Fourier transformation. 
Here k max denotes the upper cutoff of the distribution of decay constants. In any event, one 
expects a very strong decay in agreement with the behavior in w(t) for times 20/s < t < lps. 
(iii) The more gradual decrease of cr{v) for v < 10 11 is also recovered in our simulated w(t) 
curves for the two lower temperatures. For T > 950 K w(t) is roughly constant for t > lps. For 
the sodium silicate system the dispersion of a(u) as roughly expressed by o{y w 10 u Hz)/adc is 
about 10 7 at room temperature. It becomes smaller with increasing temperature and vanishes 
at around T w 900 K. 

The postulation of a crossover from vibrational-type to long-range dynamics around v = 10 11 
Hz, as derived from conductivity spectra, is consistent with physical intuition. There is, how- 
ever, no strict derivation. Here simulations may give additional information because the micro- 
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scopic nature of dynamics is known in detail. In particular, we have identified the contribution 
w local {t) which arises from local vibrational dynamics, see Fig. 12. As expected its value is close 
to w(t) for short times. The long-range contribution to w(t) can be estimated as w(t) —wi oca i(t). 
This is included in Fig. 13 for the two lowest temperatures. For these temperatures one can 
clearly see that the dispersion for t > lps is mainly due to long-range motion. For T = 640 K 
the dispersion, i.e. w(t = lps)/w(t — > oo), is around 20. 



4 Discussion and Summary 

The analysis as presented above allows us to improve our understanding of the frequency 
dependence of the conductivity. Here it is helpful that the time-dependence of w(t) (derivative 
of the mean square displacement) is strongly related to the frequency dependence of cr(v). 
As outlined above several findings like the activation energy or the Haven ratio are in good 
agreement with the corresponding results of experiments. In particular it turns out that only for 
temperatures below 900 K the system displays dispersive dynamics due to long-range dynamics. 
Thus only simulations in the nanosecond regime can cope with the effect of back- and forthjumps 
as observed in the experiment. At first view this result disagrees with the fact that even at 
T = 1240 K the mean square displacement (see Fig.4) is not diffusive but shows a subdiffusive 
regime. Our analysis of w(t), however, has revealed that the subdiffusive regime for this high 
temperature is related to local vibrational dynamics. Therefore it only contributes to the 
high-frequency regime of cr(i>), reflecting vibrational properties. 

Long-range dynamics prevails for t > 1 ps or, correspondingly, v < 10 11 Hz. This observa- 



tion seems to contradict the conclusion in Ref.||45| that the nearly constant loss and thus the 
dynamics for frequencies much lower than 10 11 Hz is due to anharmonic vibrational dynamics 
rather than jump dynamics. A priori, however, it is not clear whether long-range dynamics 
automatically implies jump-like motion and further work is necessary to clarify this point. 
As shown in Fig. 13 the subtraction of wi oca i(t) from w(t) yields a plateau-like region for short 
times. This is in agreement with recent work on amorphous 0.5 Ag2S - 0.5 GeS2 ||. Of course, 
for experimental data it is not possible to determine the vibrational contribution individually. 
For some materials, however, it can be estimated reliably M. The interpretation of w(t) — 
wiocaiit) as the contribution of jumps between discrete sites is definitely justified for long times. 
In general it is justified as long as the vibrational and the jump dynamics are uncorrelated 
to each other. Intuitively one would expect that this is the case for time scales which are 
significantly longer than the time it takes for the ion to transfer between two adjacent sites. Two 
important questions remain. What is this crossing time? What is the physical interpretation 
of w(t) — Wiocaiit) for shorter times than the crossing time? Work along this line is in progress. 
Another important aspect is related to the time-temperature superposition principle. The 
scaling, observed in Fig. 6, is consistent with the Summerfield scaling as can be easily checked 
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with the help of Eq.[|. As is evident from Eq.|I] measurement of cr(u) is not sufficient to obtain 
an absolute value of the mean square displacement. From measurements of e'(oo), however, it 
becomes possible to get information about the absolute value c(T) of (r 2 (t)) for short times 
(or high frequencies, respectively) where only the local fluctuations contribute [H, |], In 
agreement with our results (see Fig. 11) one finds to a good approximation c(T) oc T [}i9 |. A 
more detailed comparison of the experimental scaling results and our findings for the mean 
square displacement is not possible so far since our scaling regime is still very limited and 
simulations at still lower temperatures had to be taken into account. 

Finally, we would briefly comment on the different activation energies of the conductivity and 
the incoherent scattering function. In glass-forming systems at low temperatures a similar effect 
is seen when comparing the rotational correlation function (which basically behaves like the 



translational incoherent scattering function) and the diffusion constant. As explained in Ref. [fH 



such an effect may be related to a broadening of the waiting time distribution for decreasing 
temperature. The diffusion constant is mainly dominated by the fast ions. In contrast, decay 
of the incoherent scattering function (with time scale r) requires that also the slow ions move. 
Thus an increase of heterogeneity would increase the product of the diffusion constant and 
the relaxation time r. On a qualitative level this is seen in our simulations. Although the 
heterogeneities indeed increase with decreasing temperature ||7| we do not know whether this 
is the only explanation for the different activations energies. In any event, this question is 
accessible experimentally since via multidimensional NMR observables can be measured which 
are believed to behave like the incoherent scattering function ||48|| . 

The increase of heterogeneity as indirectly seen by the different activation energies implies that 
in a very strict sense the dynamics at lower temperatures is not a simple time-scaled version 
of the dynamics at higher temperatures. Thus we have to conclude that the mean square 
displacement and the incoherent scattering function, which both can be interpreted as reduced 
representations of the dynamics, are not sensitive to these deviations from time-temperature 
superposition. Actually, a similar observation has been already reported for a glass-forming 
hard sphere system for which the non-gaussian parameter strongly depends on density although 
the mean square displacement displays time-density superposition pl| . 

Several important questions remain to be answered. How do the observed features correlate 
with the local network structure (see Refs. fT8|, |T9[ | for interesting features observed recently)? 
To which degree do the observed features change with concentration? What can one learn 
about the nature of correlated back- and forthjumps? For the latter question it turns out to 
be very helpful to study appropriate three-time correlation functions |47| . 
Of course, the results of several experiments as well as previous simulations already contain 
part of the answers. However, since computer simulations are always dealing with temperatures 
much higher than experiments one has to be very careful to extract information from simulations 
which is relevant for experiments. Therefore the advent of faster computers and thus the ability 
to simulate at sufficiently low temperatures will render computer simulations of ion conductors 
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a fruitful field for the future. 
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Fig. 1 Trajectories of three different lithium ions at T = 750 K as projected on a plane. 
The box length is 9 A. 

Fig. 2 The partial structure factor g(r) for the lithium ions at different temperatures. 

Fig. 3 The integrated partial structure factor h(r) for the lithium ions at different temper- 
atures. 

Fig. 4 The mean square displacement (r 2 (t)) for lithium at different temperatures. 

Fig. 5 The temperature dependence of the diffusion constant of Li 2 SiC>2. The data 
marked as circles have been taken from Ref . [ 26[] . 

Fig. 6 The same as in Fig.4. The individual times have been scaled in order to show the 
time-temperature superposition. 

Fig. 7 The time dependence of the incoherent scattering function S(q, t) for q = q max = 
27r /do at different temperatures. The broken lines correspond to a fit with the KWW 
function y4exp(-(t/r) /3 ) with (3 = 0.45. 

Fig. 8 The time dependence of dynamic cooperativity as expressed by N coop (t) for different 
temperatures. 

Fig. 9 The self part of the van Hove correlation function G s (r,t) for different times t at 
T = 750 K. 

Fig. 10 The local mean square displacement (r 2 (t))i oca i which averages over all displacements 
less than r min = 1.5 A. 

Fig. 11 The temperature dependence of (r 2 (t — ► oo)) i oca i- 

Fig. 12 Wiocaiit) at different temperatures. 

Fig. 13 w(t) at different temperatures. For the two lowest temperatures also w(t) —wi oca i(t) 
is displayed. 
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